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CALCULATION OF IRROTATIONAL WIND PATTERN WITH 
APPLICATION TO CLEVELAND TOPOGRAPHY 
by Robert Siegel 
Lewis Research Center 

SUMMARY 

Small perturbation theory is applied to compute the deflection of the wind blowing 
across land that has an irregular topography. As an illustration, the method is applied 
first to the flow around a single hill of Gaussian profile. Then calculations are made for 
the irregular topography on the east side of Cleveland where the elevation changes by sev- 
eral hundred feet. It was found that the topography produced small wind deflections that 
would not be of practical importance in air pollution dispersion studies. The calculations 
were for a neutrally stable atmosphere. Although they are not investigated here, other 
factors such as thermal stratification of the atmosphere, diurnal variations, and convec- 
tion currents resulting from the proximity of Lake Erie and the city heat island effect are 
expected to be more significant than the influence of topography. 


INTRODUCTION 

In order to use the diffusion equation or various dispersion models to calculate the 
dispersion of pollutants from emission sources in a city, it is necessary to know the local 
wind velocities. For example in reference 1 the diffusion equation was solved numer- 
ically for the Los Angeles area. The wind velocities were interpolated from measure- 
ments at 32 stations throughout the city. In the simplified dispersion model devised by 
Hanna (ref. 2), which is applicable for distributed area sources of pollutants, the local 
concentration was found to depend inversely on the local wind velocity. Even though the 
analytical result in reference 2 has a very simple form, it was found that good concentra- 
tion results were obtained in calculations performed for the Los Angeles and Chicago 
areas. 

One of the factors influencing the local wind velocities and flow pattern is ground 
topography. In reference 3 computations were made of the smoke dispersion from an 


electric power plant. Downwind in the direction of the prevailing wind were two hills 
that influenced the wind pattern. The situation in greater Cleveland is somewhat sim- 
ilar. On the east side of the city there is a rise in the land in the direction of the wind, 
which is predominantly from the west and southwest. In this region the land also rises 
to the south away from the lake as shown by the simplified contour map in figure 1. This 
figure also shows the location of the Cuyahoga river valley in the central portion of the 
city where heavy industry is concentrated. 

The object of this report is to determine whether topography will significantly deflect 
the wind that is passing over the industrial region and then blowing across the east side 
of the city. This would have some bearing on the interpretation of information from pol- 
lution measuring stations distributed throughout the city and would provide information 
for use in dispersion computations. 

To obtain a first approximation as to whether the irregular topography would signif- 
icantly deflect the wind, it was felt that boundary -layer effects could be neglected and an 



Figure 1. - Simplified contour map of greater Cleveland area showing the increase in elevation in the southeast area. (Contour elevation lines 
are in feet, t 
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inviscid irrotational analysis performed for a neutrally stable atmosphere. If the invis- 
cid results show appreciable effects on the flow, then these results can be used as the 
outer flow and the boundary layer effects can be calculated. An inviscid calculation 
would not be valid in a region where there is flow separation such as on the downstream 
side of a hill. However, as discussed in reference 4 (p. 29), inviscid irrotational flow 
theory provides a reasonable quantitative description of the flow pattern on the windward 
side of a hill. The prevailing wind direction in Cleveland is directed toward the rising 
land in the region of interest here, and hence the inviscid analysis should be reasonable 
for a first approximation. 


SYMBOLS 

A height of undisturbed streamline above x-y plane 

H dimensionless height of ground above x-y plane, h /l r 

h height of ground above x-y plane 

h max maximum height of hill 
Z r reference length 

p pressure 

t intermediate variable, £ + It] 

u, v,w velocities in positive x, y, z directions 
u^ undisturbed wind velocity in positive x direction 
W complex potential function, 4> + i^ 

X, Y, Z dimensionless rectangular coordinates, X is parallel to direction of undisturbed 
wind 

X',Y' dimensionless grid coordinates for Cleveland 

x,y,z rectangular coordinates, x is in direction of undisturbed wind 

z complex variable, x + iz 

a angle with respect to -x axis of approaching undisturbed wind 

£ dummy z variable 

r\ dummy y variable; imaginary part of t 

£ dummy x variable; real part of t 

7} dimensionless variables, %/l r , r]/l T 

p density 
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a constant in Gaussian distribution 
potential, real part of W 
(p perturbation potential 
\f r streamfunction, imaginary part of W 

4^ streamfunction for n tn streamline 

ANALYSIS 

To gain some appreciation of how the flow is affected by topographical features, the 
flow across two configurations of regular geometry will be briefly examined: a two- 
dimensional step and a hill with a Gaussian profile. All of the calculations are for a neu- 
trally stable atmosphere and flow separation effects are neglected. 


Deflection of Flow by a Two-Dimensional Step 

The geometry for the two-dimensional step is shown in figure 2. At a great distance 

before reaching the step, the undisturbed flow velocity has components u = -u m and 

1 °° 

v =v oo so that the velocity vector is at angle a = tan' (v^-u^) relative to the negative 
x direction. To show how the flow is deflected, the streamlines will be obtained for this 
approaching flow. For an inviscid irrotational calculation it is shown in appendix A that 
the velocity in the y -direction will be v m everywhere, and the component of the velocity 
vector in the x-z plane can be treated separately from that in the y -direction. In refer- 
ence 5, the solution is given in the x-z plane for a two-dimensional flow approaching in 
the -x direction and then passing over a step. Since the effect of the v -velocity compo- 
nent is uncoupled , this solution gives the streamlines projected in the x-z plane for the 
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present case. The interest now is in obtaining a streamline path projection in the x-y 
plane by including the effect of the v-velocity component. 

Along a streamline 


dy = v 
dx u 

Since v is everywhere equal to v^ this equation can be integrated to obtain 

y( x ) = v oo / -dx (1) 

•'streamline u 

where the integration is along a streamline. To obtain u along a streamline the solution 
in reference 5 is used. This gives the complex potential W as a function of x and z 
in terms of an intermediate variable t, 


W = $ + uk = — — cosh t (2) 

v 

z=x + iz = -(t + sinh t) (3) 

v 

where h is the height of the step. Then the complex conjugate velocity in the x-z plane 
is given by 


u-iw = -^=-^^L = u sinh t 
dz dt dz 00 1 + cosh t 


(4) 


If we let t = £ + irj, taking the imaginary part of equation (2) gives 

-U oo h 

= sinh £ sin 77 (5) 

7T 

IV 

Along the n in streamline, which originates at a distance A above the base plane in fig- 
ure 2, is a constant. If we let \k n = 0 along z = 0, then \k n = -u w A since the dif- 
ference between stream function values for streamlines is equal to the volume flow be- 
tween them. Hence the relation between £ and 77 along a streamline is, from equa- 
tion (5), 
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= sinh | sin 77 
h 


( 6 ) 


To obtain u, take the real part of equation (4) to give 


u 


= 


sinh g 

cosh £ + cos 7? 


(7) 


Then u along a streamline is given in terms of 77 by eliminating £ in equation (7) by 
use of equation (6). This yields 


u 

u_ 


* 


7tA 

h 


n 



( 8 ) 


To carry out the integration in equation (1), the relation between x and 77 along the 
streamline must be known. In addition, to be able to plot the streamlines a relation be- 
tween z and 77 along a streamline is needed. These relations are obtained by taking the 
real and imaginary parts of equation (3) which gives, 


- = - [ £ + sinh | cos 77] (9) 

h 7T 


— = i [77 + cosh £ sin 77] (10) 

h jt 


Then eliminate £ by use of equation (6) to obtain x and z along a streamline in terms 
of 77 



By evaluating equations (11) and (12) for various 77, the streamlines projected in the 
x-z plane are plotted in figure 3(a). 
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The integration in x is carried out along a streamline using equation (11) to obtain 77 (x) 
so that the integrand can be evaluated at any x value. The Xq in the lower limit of in- 
tegration is a value taken far enough before the step so that the flow has not yet been de - 
fleeted (xq « 6h was reasonable in most instances). 

It is evident from equations (13) and (11) that the quantity (y/h) (-u^/vj depends 
only on the single parameter A/h which is the height of the undisturbed streamline 
above the x-y plane divided by the height of the step. The quantity (y/h) (-u^/vj is plot- 
ted as a function of x/h in figure 3(b) for streamlines with three different A/h. If 
-u^ = v^, that is, if the wind approaches the step at a = 45°, these curves are the actual 
streamline paths and show the horizontal wind deflection. For v M + -u M the horizontal 
deflection is found by multiplying the ordinate of figure 3(b) by v^-u^. 


Analysis for a Three-Dimensional Topography 

For a three-dimensional variation in land elevation, the solution for the Laplace 
equation governing the velocity potential would usually not be possible analytically and 
could be obtained numerically by finite differences or by a technique such as in refer- 
ence 6. In many instances, however, the changes in elevation are rather gradual (sev- 
eral hundred feet per mile) and consequently a small perturbation theory can be applied 
to obtain deviations in velocity from the free -stream value. Some of the basic aspects of 
the small perturbation theory are given in reference 7. This theory was applied to a real 
topography in reference 3 although the details of the analysis are not given there. The 
small perturbation equations used here are derived in appendix B. The results for the 
perturbation velocities are 


U (X,Y,Z) = 1 + 


u 


i/7 


-2(X - l) 2 + (Y - ^) 2 + Z 2 
r |=-=o [(X - h 2 + (Y - ri) 2 + Z 2 ] 5/2 


H(|,»7)d? dr? 



H(i,^)d | d^ 


y (i4) 


H(£,77)d£ dr? 
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The H(X,Y) is the local elevation of the land above the X-Y plane normalized by a char- 
acteristic dimension. In the integrals the X, Y, and Z are the coordinates of the location 
at which the velocity perturbations are being computed. The integration over £ and 7 ) 
takes into account the effect that the surrounding land contour has on the velocity 
perturbations. 

Deflection of Flow by a Hill 

As an illustration of the flow for a three-dimensional topography, the wind pattern 
around a hill will be considered. For a three-dimensional topography such as a hill, 
there is a relief effect (ref. 7) resulting from some of the flow being diverted around the 
hill rather than having to flow over the top of it as in a two-dimensional situation. A sim- 
ple hill geometry can be specified by using a Gaussian contour of revolution (all lengths 
are nondimensionalized by the maximum hill height, h max ) 

Ifey) = H(X, Y) = e ~(X 2 +Y 2 )/a 2 ( 15 ) 

max 

(Another function used in discussing wind flows (ref. 8) is the "mountain" function 
[l + (X 2 + yV/ 3 2 ]- 1 where /3 is a constant.) The or in equation (15) determines the 
steepness of the hill. Results were carried out here for a - 2 and the hill contour for 
this value is shown in figure 4. The H(X, Y) was inserted into equations (14), and the 
integrations carried out numerically for various (X, Y, Z). Some results are tabulated in 
table 1(a) for the first X,Y quadrant. This is sufficient to give the entire flow pattern 
because the inviscid flow is symmetric around the hill as flow separation effects are not 
taken into account. The velocity vectors at various heights in the X-Z plane are shown 
in figure 4(a). 

Since measurements of the wind are often at a fixed height above the ground, a set of 
calculations of this type was also performed where the height above the ground was 0. 1 
of the maximum hill height. For example, for a hill 500 feet high (typical of Cleveland 
topography) the results show the wind 50 feet above the ground. For this calculation, 
the Z(X,Y) of equation (14) was set equal to 

Z(X, Y) = H(X, Y) + height above ground = H ( X> y) + 0. 1 

''max 

Figure 4(b) shows the horizontal velocity vectors at various X,Y values and illustrates 
the flow deflection around the hill. 
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4.0 3.2 2.4 1.6 .8 0 .8 


R- X 2 + Y 2 

(b) Horizontal velocity vectors for flow around hill at fixed distance of one-tenth of maximum hill height above 
the ground Z(X, Y) - H(R) + 0. 1. 

Figure 4. - Wind velocity vectors for flow over Gaussian hill (o ■ 2). 


TABLE I. - VELOCITY COMPONENTS 
[Each set of three values is u/u^ v/u^ 

(a) Velocity components for flow over hill with Gaussian profile (a = 2) 


Dimensionless rectangular coordinate 
in direction of undisturbed wind, 

X 


0.5 0 Inside hill Inside hill 0.983 0.914 

0 0 
-.229 -.125 

1 Inside hill Inside hill 0.989 0.929 

-.068 -.046 


-.186 -.103 


1.001 

0.959 

-.083 

-.059 

-. 102 

-.058 


1.009 0.984 
-.059 -.046 

-.039 -.024 


1.007 0.955 

0 0 
-.146 -.095 


1.008 0.962 
-.039 -.030 

-.123 -.081 


67 

1.011 

0.978 

39 

-.052 

-.041 

61 

-.075 

-.051 

42 

1.012 

0.992 

30 

-.041 

-.035 

27 

-.034 

-.025 


Dimensionless rectangular coordinate 
in direction of undisturbed wind, 

X 


1.013 0.996 

-.028 -.026 
-.028 -.023 


1.016 0.990 

0 0 
-.063 -.052 


1.016 0.992 
-.014 -.013 

-.056 -.047 


37 

1.014 

0.996 

14 

-.021 

-.019 

28 

-.039 

-.034 


1.012 1.000 
-.019 -.019 

-.022 -.020 









































TABLE I. - Concluded. VELOCITY COMPONENTS 


[Each set of three values is u/u^, v/u^, w/u^.] 


(b) Velocity components at height above ground 
equal to 0 . 1 of the maximum hill height for 
flow over a Gaussian hill (a = 2) 


Y 

Dimensionless rectangular coordinate in 
direction of undisturbed wind, 

X 


0 

1 

2 

3 

4 

5 

0 

1. 145 

1. 114 

0.980 

0.887 

0.909 

0.955 


0 

0 

0 

0 

0 

0 


0 

-. 140 

-.236 

-. 146 

-.042 

-.007 

1 

1.156 

1. 115 

0.982 

0.897 

0.923 

0.960 


0 

-.044 

-.077 

-.060 

-.027 

-.011 


0 

-. 139 

-.205 

-. 117 

-.033 

-.006 

2 

1. 147 

1.097 

0.991 

0.940 

0.952 

0.972 


0 

-.077 

-.108 

-.076 

-.037 

-.017 


0 

-.102 

119 

-.058 

-.016 

-.003 

3 

1.088 

1.058 

1.004 

0.978 

0.978 

0.984 


0 

-.061 

-.077 

-.056 

-.031 

-.016 


0 

-.039 

-.038 

-.018 

-.005 

-.001 

4 

1.041 

1.029 

1.009 

0.996 

0.992 

0.992 


0 

-.029 

-.038 

-.031 

-.021 

-.013 


0 

-.008 

-.008 



-. 004 

-.001 

-.001 


Application to the Topography of Greater Cleveland 

As illustrated by the results for the Gaussian hill, the small perturbation method is 
adequate to provide flow deflections for gradual changes in land elevation. To insert the 
land contour of Cleveland into equation (14), it was decided to use a 1 mile square grid 
alined east -west and north-south and interpolate intermediate points for purposes of the 
double integrations. This grid has been designed as X',Y’ to distinguish it from the X 
and Y directions which correspond to the horizontal directions along and normal to the 
undisturbed wind. The origin of the grid was placed at downtown Public Square. The 
grid and ground elevations are shown in figure 5. Lake Erie is the region of constant 
elevation (584 ft) in the upper part of the figure. To normalize the distances, a char- 
acteristic length of 528 feet (0. 1 mile) was selected; thus each grid spacing was 10 units. 
The normalized elevation H(X, Y) above lake level, which was taken as the base plane, 
was then given by 

H(X Y) = eleva *i° n at (X, Y) - 584 
528 

The calculations for the wind were made at a level 50 feet above the ground so that the Z 
in the integrals of equation (14) is given by 
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Locations where wind was calculated 



Figure 5. - Land elevations in feet for 1 mile grid of greater Cleveland area (X 1 and Y' are in units of 1/10 mile). 











Z = H(X, Y) + J^L 
528 


In the analysis, the undisturbed flow is in the X direction. This would correspond 
to a wind from the west when X is alined with the X’ on the grid of figure 5. The pre- 
vailing wind in Cleveland is from the west or southwest. To obtain results for a south- 
west wind it is necessary to use a rotation of coordinates to obtain the ground elevations 
H(X,Y) used in the integration. For a southwest wind the X-Y integration coordinates 
are rotated counterclockwise 45° from the x'-Y* coordinates in figure 5. Hence, to ob- 
tain H(X,Y), the corresponding X', Y' are found from 

X’ = (X - Y) Y' = (X + Y) 

2 2 

and then the H is obtained at these X*, Y* in figure 5. 

The main topographical features evident on figures 1 and 5 is that the west side of 
greater Cleveland is quite flat. Relative to the surrounding land there is a decrease of 
elevation of a few hundred feet in the Cuyahoga river valley where heavy industry is lo- 
cated. The main feature is a rise in elevation of several hundred feet on the east side 
into the ’’Heights" area (Cleveland Heights and Shaker Heights). The land rises from 
west to east and from north to south so that an elevated plateau exists in the form of a 
rounded corner when viewed from above. With regard to the dispersion of pollutants, it 
is desired to know whether the wind from the west or southwest that passes over the in- 
dustrial section is significantly deflected when passing into the heights area. 

The wind velocity perturbation integrals in equation (14) were carried out at loca- 
tions X* = 40, 60, and 80 for Y* = -40, -20, 0, and 20. As shown in figure 5, these 
locations are in the region of the largest elevation change and hence provide the largest 
wind deflections. The results for the velocity perturbations are summarized in table II 
and will be discussed in the next section. It is evident, however, that the velocity per- 
turbations are quite small. 


DISCUSSION 

To obtain an indication of the flow deflection by topographical features, the stream- 
line deflection was computed for an initially uniform horizontal wind incident at an angle 
across a step as shown in figures 2 and 3 The u component of velocity decreases as 
the step is approached as shown by the spreading of the streamlines projected in the x-z 
plane. The decrease in u component provides more time for the v component to move 
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the flow along the y -direction. Consequently, during the approach to the step, the flow 
is turned along the y -direction. At the crest of the step the opposite effect occurs; the 
u component increases, and the flow is turned toward its original direction. The ordi- 
nate in figure 3(b) can be interpreted directly as the streamline deflection for an incidence 
angle of a = 45° (v w - -uj. Streamlines for any other angle can be found by multiplying 
the ordinate in figure 3(b) by v^-u^ = tan a. The parameter A/h is the fraction of the 
step height at which the streamline of the incoming undisturbed flow is located above the 
x-y base plane. The streamlines near the ground are deflected the most as a result of 
having the largest decrease in u velocity when approaching the step. The results show 
that the deflections for incidence at 45° are of the order of a few step heights; when 
passing over a 500 -foot cliff, the wind originating 50 feet from the ground would be hor- 
izontally displaced about 1000 feet. 

For a three-dimensional flow around a hill there is a relief effect with regard to the 
acceleration of flow over the crest of the hill. This is because the flow is partially di- 
verted around the sides of the hill. The flow is also decelerated somewhat as it ap- 
proaches the hill in a fashion similar to the approach to a step. These two factors con- 
tribute to the decrease in the u velocity (the component in the incident flow direction) 
for the flow approaching the hill and produce the values of u/u M that are less than 1 in 
table I. To illustrate the three-dimensional calculation technique, a few results were 
carried out for a hill of Gaussian cross section (fig. 4). The results in figure 4(b) are 
for a height above the ground of one -tenth of the maximum hill height. The turning of the 
flow is always small; the perturbations of the velocity components (table I) are in the 
range of several percent. 

Since the variations in greater Cleveland land elevation are in the range of several 
hundred feet, the preceding calculations for simple geometries would lead one to expect 
that the deflections of the wind by the topography would be rather small. Two sets of 
calculations were carried out, and both substantiate this expectation. The velocity per- 
turbations for flow across the east side of the city were at most only a few percent of the 
mean velocity in the regions of greatest change in ground elevation. This was true for 
winds from the west and from the southwest. The w/u^ components are approximately 
equal to the slope of the land in the x- direction as would be expected for gradual changes 
in contour. Note that u is in the direction of the undeflected wind and v and w are 
normal to that direction. Since the velocity perturbations are small, it is felt that the 
topographical features would represent a second-order effect on the wind deflection as 
compared with other factors. Factors of possible significance are convection currents 
resulting from temperature inequalities arising from the lake adjacent to the city, heat 
island effects, diurnal temperature variations, and thermal stratification of the atmos- 
phere. These would provide complicated transient wind variations. 
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CONCLUSIONS 


A small perturbation technique was applied to compute the three-dimensional veloc- 

> 

ity perturbations for wind blowing across the topography of greater Cleveland. The cal- 
culations were irrotational and for a neutrally stable atmosphere. The velocity perturba- 
tions were found to be small, indicating that the small perturbation technique used in the 
analysis was adequate. It is concluded that the topographical features of greater Cleve- 
land provide only a second-order effect on the wind as compared with other factors of 
possible significance such as diurnal temperature variations and heat island effects. 

Lewis Research Center, 

National Aeronautics and Space Administration 
Cleveland, Ohio, August 21, 1972, 

501-24. 
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APPENDIX A 


INViSCID IRR0TAT10NAL EQUATIONS 


For inviscid flow in three dimensions, the Euler equations are 


ulH + vlH+w^i= -IIP 

3x 3y 3z p 3x 


uiI + vlI+w9l=-i®E 

3x 3y 3z p 3y 


u^! + v— + w— = -1^. 

3x 3y 3z p 3z J 


(Al) 


For flow across the two-dimensional step in figure 2, there is no variation in the y- 
direction so that 3/3y =0. Then equations (Al) become 


u — +w— = -I& 
3x 3z p 3x 


3x 3z 


> 


u 3w + w Sw = _iap 
3x 3z p 3z J 


(A2) 


For irrotational motion 


3v_3u_0 3w _ 3v _ o 3u 3w _ q 

3x 3y 3y 3z 3z 3x 

and with 3/3y = 0 these yield 3v/3x = 0, and 3v/3z = 0. Hence, v must be a constant 
throughout the flow and is v^ everywhere. The Euler equations (eqs. (A2)) then further 
reduce to 


18 



(A3) 


u + w — = - > 

3x 3z p 3x 

u 3w + w 3w = _ 1 3p 
3x 3z p 3z 

which are the two-dimensional equations for u and w in the x-z plane. The u and w 
components are thus uncoupled from the constant v=v w component. 
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APPENDIX B 


SMALL PERTURBATION ANALYSIS FOR THREE- 
DIMENSIONAL VELOCITY DISTRIBUTION 

The quantities u, v, and w are the velocity components, and for the situation con- 
sidered here ju -u^l « ju^J, |v| « |uj, and |w| « [u^. Let the velocity poten- 
tial be $(x,y,z) = u M x + <p(x,y, z) where cp is the perturbation potential. Then 3<j>/3x = 
u oo + 3^/3 x = u, 33>/3y = dcp/dy = v, and 3<J>/3z = 3<p/3z = w. The local elevation of the 
ground above the x-y base plane is given by h(x,y) (fig. 6). At the surface of the ground 
the boundary condition requires that the streamlines are tangent to the surface. Hence 
at the surface 


u _ v _ w 
dx dy dz 

Then within the approximation of small perturbation theory 


dz r / v t dz dz 
W = U — = Ko + (u - Ujj — = Uoo — 
dx dx dx 


surface 


Using w = 3<p/3z and the fact that dz/dxj surface = 3h/3x gives 


= u — 
3z °° 3x 


(Bl) 


Equation (Bl) expresses the derivative normal to the x-y plane of the unknown pertur- 
bation potential in terms of the known free -stream velocity and slope in the undisturbed 
flow direction of the land contour . For small perturbation theory this boundary condition 
at the surface of the ground is applied at the base plane z = 0 (fig. 6). 

The solution for the perturbation potential requires solving Laplace's equation 
o 

V (p = 0 in the half space 0 < z < °° subject to the boundary condition that the normal 
derivative at z = 0 is given by equation (Bl). The Greens function for this problem is 
given in reference 9 as 
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g = 


lim 

C—0* 


[(I - x) 2 + (n - y) 2 + (K - z) 2 ] 1/2 [u - x) 2 + ( V - y) 2 + (c + z) 2 ] l/2 


[(I - x) 2 + (?? - y) 2 + z 2 ] 1 / 2 
Then from Greens fundamental solution (ref. 9) 


(B2) 


<p(x 


,y, z) = - — f f g-^i?dCd77 
47 r 4 / 3? 


x-y 

plane 


Substituting equations (Bl) and (B2) gives 


<Kx,y,z) 



1 

[it -X ) 2 + (77 - y) 2 + z 2 ] 1//2 


— d£ dr) 


(B3) 
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Integrate by parts with respect to i; to obtain 


[(|-x) 2 + (7 7 -y) 2 + z 2 ] 1/2 | 


/ 


- i)2(| - x)h 


d - x) 2 + (7? - y) 2 + z 2 ] 3 / 2 


d^>d?7 


cp(x,y,z) = 


n co 


(x - g)h(g,i7) 


[(x - £) 2 + (y - ry) 2 + z 2 ] 3 / 2 


d£ d77 


Divide all lengths by a characteristic dimension l r to obtain the dimensionless form 


$X,Y,Z) = 


n c 

» 


g - IsM di ai ( B4) 

(X - £) 2 + (Y - J)) 2 + Z 2 ] 3/2 


The velocity components are given by 


u(x,y,z) _ 1 + J_ 8<p(x,y,z) = 1 + dtp 
u„ 8x 3X 


v(x,y, z) _ J_ d<p(x, y, z) _ 8^> 
u =o u oc 3 y 9Y 


w(x,y,z) _ 1 d<p(x,y,z) _ dtp 


u u 8z 

oo oo 
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By differentiating cp the velocity components are obtained as 


1 , i/; _ 

1 + _ 

U oc 

2rr 

v(X,Y,Z) _ 

_ 

U oo 

2rr .> 

w(X,Y,Z) _ 

- 

U C 0 

2 ir 


n co 


-2(X - j) 2 + (Y - y) 2 + Z 2 
(X - £) 2 + (Y -^) 2 + z 2 ] 5/2 


p OO /lOO 

v_oo V-OO 


A A. 


(X - 4)(Y - r>) 


(X - £) 2 + (Y - ^) 2 + Z 2 ] 5/2 


poo pOO 

t/_O0 «/, -oo 


X-J 

(X - i) 2 + (Y - i) 2 + Z 2 ] 5/2 


H(i, i?)di d77 


H(i,5)did5 > 


H(i,i)di dr? 
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